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I. INTRODUCTION 

Within the standard model (SM) Q, spontaneous 
breaking of electroweak symmetry gives mass to the W 
and Z bosons 0], and to the fundamental fermions via 
their Yukawa interactions with the Higgs field. In the 
SM, the symmetry-breaking mechanism predicts the ex- 
istence of one neutral scalar particle, the Higgs boson, 
whose mass (mjy) is a free parameter. 

Precision electroweak data, including the recently up- 
dated measurements of the VF-boson and top-quark 
masses from the CDF and DO Collaborations [3, J . yield 
an indirect constraint on the allowed mass of the Higgs 
boson, niH < 152 GeV/c 2 Q, at the 95% confidence level 
(C.L.) @. Direct searches at LEP2 exclude SM Higgs 
boson masses below 114.4 GeV/c 2 0. The ATLAS and 
CMS Collaborations at the Large Hadron Collider (LHC) 
have recently reported the observation of a new boson 
with mass of around 125 GeV/c 2 @, 0. Much of the 
sensitivity of the LHC searches comes from gluon-gluon 
fusion (gg — > H) production and Higgs boson decays to 
77, ZZ, and W + W~ . Published searches for associated 
production VH -t Vbb at the LHC, where V = W or 
Z pjj [H| , have not yet reached sensitivity to SM Higgs 
boson production. The CDF and DO Collaborations have 
recently reported combined evidence for a particle, with 
a mass consistent with that of the new boson observed 
at LHC, produced in association with a W or Z boson 
and decaying to a bottom-antibottom quark pair |l2| ■ 

In this article, we combine the most recent results 
of SM Higgs boson searches in pp collisions at y/s = 
1.96 TcV using the full Tevatron Run II integrated lu- 
minosity of up to 10 fb _1 per experiment. The analyses 
combined here seek signals of Higgs bosons in the mass 
range 90-200 GeV/c 2 , produced in association with a 
vector boson (qq — > VH), in association with top quarks, 
through gluon-gluon fusion, and through vector boson fu- 
sion (VBF) (qq q'q'H). The Higgs boson decay modes 
studied arc H -> bb, H -> W+W~, H -> ZZ, H -> 
r + r _ , and H 77. For Higgs boson masses greater 
than 130 GeV/c 2 , searches for H — > W + W~ decays with 
subsequent leptonic W decays provide the greatest sensi- 
tivity. Below 130 GeV/c 2 , sensitivity comes mainly from 
associated VH production, with the H boson decaying to 
bb and the W or Z boson decaying leptonically. While we 
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present our results in the full mass range, we also focus 
specifically on the mass hypothesis to# = 125 GeV/c 2 , 
due to the recent LHC findings. Specifically, we show the 
sensitivity of the searches over the full mass range to a 
SM Higgs boson signal with m# = 125 GeV/c 2 . Previous 
Tevatron SM combinations, focused respectively on the 
H — > bb and H — > W + W~ decay modes, are published 
in Refs. [HI, [H|. The results presented here are based on 
the combinations of the searches from each experiment 
as published in Refs. [hi. ITsT]. 

This article is structured as follows. Section [TT| dis- 
cusses the simulation methods used to predict the yields 
from the signal and SM background processes. Sec- 
tion |TlT] briefly describes the CDF and DO detectors. Sec- 
tion [IV] describes the event selections used by the various 
analyses and Section [V] presents the data. Section IVII 
provides a brief introduction to the statistical procedures 
used and Section IVIII discusses the different sources of 
systematic uncertainties and how they are controlled. 
Sections I Villi and IIXI present the results in the contexts 
of the SM and extensions to it. Section [X] summarizes 
the article. 



II. EVENT SIMULATION 

Higgs boson signal events are simulated using the 
leading-order (LO) calculation from pythia [la ], with 
CTEQ5L (CDF) and CTEQ6L1 (DO) [T^ parton dis- 
tribution functions (PDFs). The normalization of these 
Monte Carlo (MC) samples is obtained using the highest- 
order cross-section calculation available for the corre- 
sponding production process. The cross section for 
the gluon-gluon fusion process is calculated at next-to- 
next-to-leading order (NNLO) in quantum chromody- 
namics (QCD) with soft gluon rcsummation to next- 
to-next-to- leading- log (NNLL) accuracy 0, The 
WH and ZH cross-section calculations are performed 
at NNLO precision in QCD and next-to-leading-order 
(NLO) precision in the electroweak corrections [20]. The 



VBF cross section is computed at NNLO in QCD [21 



and the electroweak corrections are computed with the 
HAWK program [22||. All signal production cross sec- 
tions are computed using the MSTW2008 PDF set pi. 
The PDF uncertainties are assessed following Refs. [24l . 
I25] ]. The Higgs boson branching fractions are from 
Ref. [26[ and rely on calculations using HDECAY [27} and 
prophecy4f [28( . The distribution of the transverse mo- 
mentum (jpt) of the Higgs boson in the PYTHiA-generated 
gluon-fusion sample is reweighted to match the pt as cal- 
culated by HQT [29|, at NNLL and NNLO accuracy. 

We model SM and instrumental background processes 
using a mixture of MC and data-driven methods. In the 
CDF analyses, backgrounds from SM processes with elec- 
troweak gauge bosons or top quarks are modeled using 

PYTHIA, ALPGEN [3(|, MC@NLO [5l|, and HERWIG [H| . 

For DO, these backgrounds are modeled using PYTHIA, 
ALPGEN, and SINGLETOP [Hi. An interface to pythia 
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provides parton showering and hadronization for genera- 
tors without this functionality. 

Diboson ( WW, WZ, ZZ) MC samples are normalized 
using the NLO calculations from mcfm [34j . For top- 
quark-pair production (ti), we use a production cross 
section of 7.04 ± 0.49 pb HU, which is based on a 
top-quark mass of 173 GeV/c 2 and MSTW 2008 
PDFs [2^ . The single-top-quark production cross sec- 
tion is taken to be 3.15 ± 0.31 pb [36[. For many analy- 
ses, the V+jet processes are normalized using the NNLO 
cross section calculations of Ref. [37| , though in some 
cases data-driven techniques are used. Likewise, the 
normalization of the instrumental, multijet and, for the 
CDF searches, the V+heavy-flavor jet backgrounds [38[ 
are constrained from data samples where the expected 
signal-to-background ratio is several orders of magnitude 
smaller than in the search samples. For the DO searches, 
the V-flight-flavor is normalized to data in a control re- 
gion, and the V-fheavy-flavor normalization, relative to 
the y+light-flavor, is taken from MCFM. In addition, for 
the DO searches, prior to fe-tagging [39[ V+jets samples 
are compared to data and corrections applied to mitigate 
any discrepancies in kinematic distributions. 

All MC samples are processed through a geant [4(| 
simulation of the detector, and reconstructed in the same 
way as data. The effects of instrumental noise and ad- 
ditional pp interactions arc modeled using MC in the 
CDF analyses, while recorded data from randomly se- 
lected beam crossings with the same instantaneous lumi- 
nosity profile as data are overlaid on to the MC events 
in the DO analyses. 

For the H — s- W + W~ analyses, the dominant irre- 
ducible background process is diboson production, while 
the dominant reducible backgrounds are Z/7*+jets, ti, 
W+7, VT+jets, and multijet production where in the lat- 
ter three cases photons or jets can be misidentified as 
leptons. For the analyses targeting H — > bb the main 
backgrounds originate from y+hcavy-flavor-jets and ti 
production. 



events are selected by singlc-lcpton and dilcpton triggers, 
a gain in efficiency of up to 20%, depending on the chan- 
nel, is achieved by including events that pass lepton+jets 
and lepton+^r triggers. 

High-quality electron candidates are identified by as- 
sociating chargcd-particlc tracks with deposits of energy 
in the electromagnetic calorimeters when both measure- 
ments are available. High-quality muon candidates are 
identified by associating tracks with hits in the muon 
detectors surrounding the calorimeters in the CDF and 
DO detectors. Lepton candidates are categorized based 
on the quality of the contributing measurements. Tight 
selection requirements yield samples of leptons with 
low background rates from hadrons or jets of hadrons 
misidentified as leptons. Looser requirements are de- 
signed to increase the acceptance for lepton candidates 
with poorly measured or partially missing information, 
with resulting higher rates for backgrounds. To optimize 
the sensitivity of the combined results, events that are se- 
lected with high-quality leptons are analyzed separately 
from those with low-quality leptons. 

Jets are clustered from energy deposits in the electro- 
magnetic and hadronic calorimeters and, in some analy- 
ses, combine information from charged particle tracks to 
improve purity or energy resolution. The transverse en- 
ergy vector Et of a calorimeter energy deposit is E sin 9h, 
where E is the measured energy, 9 is the angle with re- 
spect to the proton beam axis of a line drawn from the 
collision point to the energy deposit, and h is a unit vec- 
tor in the plane perpendicular to the beam pointing along 
that line. The missing transverse energy $t is the mag- 
nitude of the vector opposite to the sum of the Et vec- 
tors measured in the calorimeter, after propagation of 
all corrections to the calorimetric objects and for identi- 
fied muons (which deposit only small amounts of energy 
in the calorimeters) contributing to the signal topology. 
Further details of the object reconstruction algorithms 
used in the Higgs boson searches can be found in the ref- 
erences for the individual analyses (see Tables U and [TTJ) . 



III. DETECTORS AND OBJECT 
RECONSTRUCTION 

The CDF and DO detectors have central trackers sur- 
rounded by hermetic calorimeters and muon detectors 
and are designed to study the products of 1.96 TeV 
proton-antiproton collisions [U |42[ . Most searches com- 
bined here use the complete Tcvatron data sample, which 
corresponds to up to 10 fb _1 depending on the experi- 
ment and the search channel, after data-quality require- 
ments. The online event selections (triggers) rely on fast 
reconstruction of combinations of high-py lepton candi- 
dates, jets, and missing transverse energy (J5t), defined 
below. To maximize sensitivity, all events satisfying any 
trigger requirement from the complete suite of triggers 
used for data taking are considered whenever possible. 
For instance, while most of the H —> W + W~ candidate 



IV. EVENT SELECTION 

Event selections are similar in the CDF and DO analy- 
ses, typically consisting of a preselection based on event 
topology and kinematics. Multivariate analysis (MVA) 
techniques [43| are used to combine several discriminat- 
ing variables into a single final discriminant that is used 
in the statistical interpretation to compute upper lim- 
its, p-values, and fitted cross sections. Each channel 
is divided into exclusive sub-channels according to var- 
ious lepton, jet multiplicity, and 6-tagging characteriza- 
tion criteria. This procedure groups events with similar 
signal-to-background ratio to optimize the overall sen- 
sitivity. Such subdivision allows, for example, the effi- 
cient use of poorly reconstructed leptons or those in the 
forward region, the exploitation of the different domi- 
nant signal and backgrounds when training the MVAs 
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separately in each sub-channel, or reduction of the im- 
pact of systematic uncertainties. The MVAs are trained 
separately at each value of run in their respective mass 
ranges, in 5 GeV/c 2 steps. 

For the analyses exploiting the H — > bb decay, b- 
tagging and dijet mass resolution are of great impor- 
tance. Both collaborations have developed multivariate 
approaches to maximize the performance of the 6-tagging 
algorithms. The CDF 6-tagging algorithm is based on an 
MVA [Jij], and depending on the chosen operating point 
provides 6-tagging efficiencies of 50%-70% with misiden- 
tification rates for light (u, d, s, and gluon) jets of 0.5%- 
6%. In the DO analyses, the MVA builds and improves 
upon the previous neural network 6-tagger p5l l4q ] and 
achieves identification efficiencies of about 80% (50%) for 
b jets for a light jet misidentification rate of about 10% 
(0.5%). 

The decay width of the SM Higgs boson is predicted to 
be much smaller than the experimental dijet mass reso- 
lution, which is typically 15% of the mean reconstructed 
mass. A SM Higgs boson signal would appear as a broad 
enhancement in the reconstructed di-6-jct mass distribu- 
tion. The CDF and DO Collaborations search for H — >• bb 
produced in association with a leptonically decaying W 
boson, or a leptonically or invisibly decaying Z boson. 
CDF also contributes searches for WH + ZH — > jjbb 
and tiH — > tibb, where in the latter case one of the top 
quarks decays to a leptonically decaying W boson. 

Both collaborations search for the H —> W + W~ sig- 
nal in which both W bosons decay leptonically by se- 
lecting events with large missing transverse energy and 
two oppositely-charged, isolated leptons. The presence 
of neutrinos in the final state prevents reconstruction of 
the Higgs boson mass. Other observables are used for 
separating the signal from background. For example, the 
azimuthal angle between the leptons in signal events is 
smaller on average than that in background events due 
to the scalar nature of the Higgs boson and parity viola- 
tion in decays. Furthermore, the missing transverse 
momentum is larger and the total transverse energy of 
the jets is lower than they are typically in background 
events. The DO Collaboration also includes channels in 
which one of the W bosons in the H — > W + W~ process 
decays leptonically and the other hadronically. 

Although the primary sensitivity at low mass (ran < 
130 GeV/c 2 ) is provided by the H — > bb analyses and 
at high mass (m H > 130 GeV/c 2 ) by the H -> W+W~ 
analyses, significant additional sensitivity is achieved by 
the inclusion of other channels. Both collaborations 
contribute analyses searching for Higgs bosons decaying 
into tri-lepton final states, tau-lepton pairs and diphoton 
pairs. The full list of channels included is shown in Ta- 
bles Hand |H] which summarize, for the CDF and DO anal- 
yses respectively, the integrated luminosities, the Higgs 
boson mass ranges over which the searches are performed, 
and references to further details for each analysis. 



V. CANDIDATE DISTRIBUTION 

The number of contributing channels is large, and sev- 
eral different kinds of discriminating varibles are used. 
Visual comparison of the observed data with the predic- 
tions is challenging in some of the sub-channels due to low 
data counts. For a more robust comparison, we display 
the data from all the sub-channels together, aggregat- 
ing bins with similar signal to background ratios (s/b) 
from all contributing sub-channels. We collect the signal 
predictions, the background predictions, and the data in 
narrow bins of s/b, summing the contributions from bins 
in the final discriminant histograms in the sub-channels. 
A fit of the background model (see Section IVI[) to the 
data is performed before this aggregation procedure, in 
order to provide the best prediction for the background 
model in bins with the highest sensitivity. The classifi- 
cation of analysis events according to their s/b preserves 
the importance of each of the events in the histogram, 
to the extent that they arc not added to other events 
that are selected with different s/b. This representation 
of the data is not used to compute the final results, since 
the distribution indiscriminately sums unrelated back- 
grounds which are fit separately. It docs, however, pro- 
vide a guide to how much individual events contribute 
to the results and how well the signal is separated from 
backgrounds in the combined search. The resulting dis- 
tribution of log 10 (s/6) is shown for run = 125 GeV/c 2 in 
Fig. [T] demonstrating agreement with background over 
five orders of magnitude. 



VI. STATISTICAL TECHNIQUES 

The results are interpreted using both Bayesian and 
modified frequentist technique s, sepa rately at each value 
of tuh, as was done previously 0E1,|47|- The two meth- 
ods yield results that are numerically consistent; limits 
on the Higgs boson production rate typically agree within 
5% at each value of run, and with a 1% deviation when 
averaged over all positive and negative departures. For 
simplicity, when summarizing the results, we quote one 
set of values as the default, and the a priori decision 
made for the earlier Tevatron combinations to use the 
Bayesian method is retained here. Both methods use the 
distributions of the final discriminants, and not only the 
total event counts passing selection requirements. 

Each of the techniques is built on a combined likeli- 
hood (including prior probability densities on systematic 
uncertainties, 7r(0)) based on the product of likelihoods 
for the individual channels, each of which is a product 
over histogram bins, 

C(R,s,b\n,0)xH0) = I[l[ £f-^T x Ii e ~ 9 ' /2 < 
»=i j=i u ' fe=i 

(1) 

where the first product is over the number of channels 
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TABLE I: Luminosities, explored mass ranges, and references for the different processes and final states (l — e or /i, and Thad 
denotes a hadronic tau-lepton decay) for the CDF analyses. The generic labels "lx", "2x", "3x", and "4x" refer to separations 
based on lepton or photon categories. The analyses are grouped in five categories, corresponding to the Higgs boson decay 
mode to which the analysis is most sensitive: H — > bb, H — > W + W~ , H — > t+t~ , H — > 77, and H — > ZZ. 



Channel 




Luminosity 


77i# range 


Reference 






(fb- 1 ) 


(GeV/c 2 ) 




W H — y £vbb 2-jct channels 4x(5 6-tag categories) 




9.45 


90-150 


[491 


WH — > ivbb 3-jct channels 3x(2 6-tag categories) 




9.45 


90-150 


[491 


ZH — > vvbb (3 6-tag categories) 




9.45 


90-150 


[501 


ZH —> l + t~bb 2-ict channels 2x(4 6-tag categories) 


H ^bb 


9.45 


90-150 


[511 


ZH — > £+£~bb 3-jet channels 2x(4 6-tag categories) 




9.45 


90-150 


[51] 


WH + ZH -> j j66 (2 6-tag categories) 




9.45 


100-150 


[52] 


tti? — > W + bW~bbb (4 jets. 5 jcts,>6 jets) x (5 6-tag categories) 




9.45 


100-150 


[53] 


H -> W+W~ 2x(0 jets)+2x(l jet)+lx(>2 jets)+lx (low-m«) 




9.7 


110-200 


[54] 


ff -> W + VK~ (e-r had )+(M-r had ) 




9.7 


130-200 


[54] 


W.ff — > WW + W~ (same-sign leptons)+(tri-leptons) 


H -> W+W- 


9.7 


110-200 


[54] 


WH -> WW+W~ (tri-leptons with 1 T had ) 




9.7 


130-200 


[54] 


Zif ZW+W~ (tri-leptons with 1 jet,>2 jets) 




9.7 


110-200 


[54] 


H -> t+t" (1 jct)+(>2 jets) 


H -> T+r" 


6.0 


100-150 


[55] 


i? -> 77 lx(0 jct)+lx(>l jct)+3x (all jets) 


iJ — >■ 77 


10.0 


100-150 


[56] 


H ^ ZZ (four lcptons) 


# -> 


9.7 


120-200 


[57] 



TABLE II: Luminosities, explored mass ranges, and references for the different processes and final states (£ = e or fi, and Thad 
denotes a hadronic tau-lepton decay) for the DO analyses. The generic labels "lx", "2x", "3x", and "4x" refer to separations 
based on lepton, photon or background characterization categories. The analyses are grouped in four categories, corresponding 
to the Higgs boson decay mode to which the analysis is most sensitive: H — > bb, H — > W + W~ , H — > t+t~ , and H — > 77. 



Channel 




Luminosity 
(fb" 1 ) 


mji range 
(GeV/c 2 ) 


Reference 


WH —¥ £vbb 2-jet channels 2x(4 6-tag categories) 




9.7 


90-150 


[58, 59] 


WH — > £vbb 3-jet channels 2x(4 6-tag categories) 




9.7 


90-150 


[58, 59] 


ZH —> vvbb (2 6-tag categories) 


H ^bb 


9.5 


100-150 


[46] 


ZH — !• £+£~bb 2x (2 6-tag) x (4 lepton categories) 




9.7 


90-150 


[60, 61J 


H -> W+W~ -> l+vl+v 2x (0 jets,l jet,>2 jets) 




9.7 


115-200 


[62] 


H + X -> W+W- -> H+vr^ ad v (3 r categories) 




7.3 


115-200 


[63] 


H -> W+W- -> ivjj 2x (2 6-tag categories) x (2 jets, 3 jets) 
VH -> e ± Ai ± + X 


H -> VF+VK- 


9.7 
9.7 


100-200 
100-200 


[59] 
[64] 


VH £££ + X (/z/ze, 3 x e/z/z) 




9.7 


100-200 


[64] 


VH^lv nn 2x(>4jets) 




9.7 


100-200 


[59] 


VH -> r had r had /x + X (3 t categories) 




8.6 


100-150 


m 


H+X^£+T+^ d jj 2x (3 r categories) 


9.7 


105-150 


[65] 


H 77 (4 categories) 


-ff — > 77 


9.6 


100-150 


[66] 



(Nc) and the second product is over histogram bins con- 
taining riij events, binned in ranges of the final discrim- 
inants used for the individual analyses. The predictions 
for the bin contents are /i^ = R X Sij(0) + bij{9) for 
channel i and histogram bin j, where Sjj and 6jj rep- 
resent the expected SM signal and background in the 
bin, and R is a scaling factor applied to the signal. By 
scaling all signal contributions by the same factor we as- 
sume that the relative contributions of the different pro- 



cesses at each itih are as predicted by the SM. System- 
atic uncertainties are parametrized by the dependence of 
Sij and b^ on 9. Each of the n sys components of 9, 9k, 
corresponds to a single independent source of systematic 
uncertainty scaled by its standard deviation, and each 
parameter may affect the predictions of several sources 
of signal and background in different channels, thus ac- 
counting for correlations. Gaussian prior densities are 
assumed for the nuisance parameters, truncated to en- 
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FIG. 1: (color online). Distribution of log 10 (s/fe), for the 
data from all contributing Higgs boson search channels from 
CDF and DO, for m H = 125 GeV/c 2 . The data are shown 
with points, and the expected signal is shown stacked on top 
of the backgrounds, which are fit to the data within their 
systematic uncertainties. The error bars shown on the data 
correspond in each bin to the square root of the observed 
data count. Underflows and overflows are collected into the 
leftmost and rightmost bins, respectively. 



sure that no prediction is negative. 

In the Bayesian calculation, we assume a uniform prior 
probability density for non-negative values of R and inte- 
grate the likelihood function multiplied by prior densities 
for the nuisance parameters to obtain the posterior den- 
sity for R. The observed 95% credibility level upper limit 
on R, i?g5 8 , is the value of R such that the integral of the 
posterior density of R from zero to Rg§ a corresponds to 
95% of the integral of R from zero to infinity. The ex- 
pected distribution of i?gs is computed in an ensemble 
of simulated experimental outcomes assuming no signal 
is present. In each simulated outcome, random values 
of the nuisance parameters are drawn from their prior 
densities. A combined measurement of the cross section 
for Higgs boson production times the branching fraction 
B(H — > XX), in units of the SM production rate, is 
given by i? fit , which is the value of R that maximizes the 
posterior density. The 68% credibility interval, which 
corresponds to one standard deviation (s.d.), is quoted 
as the smallest interval containing 68% of the integral of 
the posterior. 

We also perform calculations with the modified fre- 
quentist technique CL S [47||, using a log-likelihood ra- 
tio (LLR) as the test statistic: 



where p(data|s + b) and p(data|6) are the probabilities 
that the data (either simulated or experimental data) 
are drawn from distributions predicted under the signal- 
plus-background and background-only hypotheses, re- 
spectively. The probabilities p are computed using the 
best-fit values of the parameters 6k, separately for each 
of the two hypotheses [67} ■ The use of these fits extends 
the procedure used at LEP [68|, improving the sensitiv- 
ity when the expected signals are small and the uncer- 
tainties on the backgrounds are large. The CL S tech- 
nique involves computing two p-values, CLb = p(LLR 
> LLR bs|fr), where LLR Q b s is the value of the test 
statistic computed for the data, and CL s+ b = p(LLR> 
LLR b s |s + 6). To compute limits, we use the ratio of 
p-valucs, CL S = CL s+ b/CLb- If CL S < 0.05 for a par- 
ticular choice of the signal-plus-background hypothesis, 
parametrized by the signal scale factor R, that hypothe- 
sis is excluded at least at the 95% C.L. The value of -Rgj? 13 
in the CL S method is the smallest value of R excluded at 
the 95% C.L. The expected limit is computed using the 
median LLR value expected in the background-only hy- 
pothesis. Systematic uncertainties are included by fluc- 
tuating around their Gaussian priors the predictions for 
Sij and bij when generating the pscudoexperimcnts used 
to compute CL s+ b and CLb- 



VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are evaluated for each final 
state, background, and signal process. Uncertainties that 
modify only the normalization and uncertainties that 
change the shape of the final discriminant distribution are 
included. To study the shape uncertainties on the distri- 
butions of the final discriminants, the relevant parameter 
is varied within one standard deviation of its uncertainty 
and the full analysis repeated using the modified distri- 
bution. For example, for the jet energy scale and reso- 
lution, the parameters of the energy scale and resolution 
are varied within one s.d. of their uncertainties and the 
analysis carried out using the kinematic distributions of 
the modified jets, also including the changes in sample 
composition resulting from the change in the jet energy 
parameters. No retraining of the MVAs is performed 
during the propagation of systematic uncertainties to the 
distributions of the discriminants. Correlations between 
signal and background, across different channels within 
an experiment and across the two experiments are taken 
into account. Full details on the treatment of the sys- 
tematic uncertainties in the individual channels can be 
found in the relevant references. 

The uncertainties on the inclusive signal production 
cross sections are estimated from the variations in the 
factorization and renormalization scale, which include 
the impact of uncalculated higher-order corrections, un- 
certainties due to PDFs, and the dependence on the 
strong coupling constant, a s , as recommended by the 
PDF4LHC working group [H| . The resulting uncertain- 
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ties on the inclusive VH and VBF production rates are 
taken to be 7% and 5%, respectively [20j]. Uncertainties 
on the branching fractions are taken from Ref . [6!| . 

For analyses focusing on gg — > H production that di- 
vide events into categories based on the number of re- 
constructed jets, the uncertainties associated with the 
renormalization and factorization scale are estimated fol- 
lowing Ref. p7pj . By propagating the uncorrelated uncer- 
tainties of the NNLL inclusive [H, [H], NLO > 1 jet 
[25j . and NLO > 2 jets [7l| cross sections to the exclusive 
gg —> H + jet. > 1 jet, and > 2 jets rates, an uncer- 
tainty matrix containing correlated and uncorrelated un- 
certainty contributions between exclusive jet categories is 
obtained. The total uncertainty on gg — > H production 
originating from these contributions varies from 10% to 
35% in individual channels depending on the number of 
jets in the final state. The PDF uncertainties are evalu- 
ated following Refs. [3[2|. 

Significant sources of uncertainty for all analyses are 
the integrated luminosities used to normalize the ex- 
pected signal yield and MC-based backgrounds, and the 
cross sections for the simulated backgrounds. For the for- 
mer, uncertainties of 6% (CDF) and 6.1% (DO) are used, 
with 4% arising from the inelastic pp cross section which 
is taken to be 100% correlated between CDF and DO. 
Cross-section uncertainties of 6% and 7% are used for di- 
boson and ti production respectively. The uncertainty on 
the expected multijet background in each channel is dom- 
inated by the statistics of the data sample from which it 
is estimated and varies from 10% to 30%. 

Sources of systematic uncertainty that affect both the 
normalization and the shape of the final discriminant dis- 
tribution include jet energy scale (1-4)%, jet energy res- 
olution (1-3)%, lepton identification, trigger efficiencies, 
and 6-tagging. Uncertainties on lepton identification and 
trigger efficiencies range from 2% to 6% and are applied 
to both the signal and MC-based background predictions. 
These uncertainties are estimated from data-based meth- 
ods separately by CDF and DO, and differ based on lep- 
ton flavor and identification category. The 6-tag efficien- 
cies and mistag rates are similarly constrained by auxil- 
iary data samples, such as inclusive jet data or ti events. 
The uncertainty on the per-jet 6-tag efficiency is approx- 
imately 4%, and the mistag uncertainties vary between 
7% and 15%. 

For the analyses targeting the H — > 66 decay, the 
largest sources of uncertainty on the dominant back- 
grounds are the rates of U+heavy flavor jets, which are 
typically 20-30% of the predicted values. Using con- 
straints from the data, the uncertainties on these rates 
are typically 8% or less. The data samples in the V+jets 
selections prior to 6-tagging are used as control samples 
to constrain systematic uncertainties in the MC model- 
ing of the energies and angles of jets. Any residual dis- 
crepancy coming from the difference between light- and 
heavy-flavor components is shown to be smaller than the 
systematic uncertainties associated with the generator or 
the correction procedures themselves. 



VIII. RESULTS - STANDARD MODEL 
INTERPRETATION 

A. Diboson Production 

To validate our background modeling and methodol- 
ogy, independent measurements of SM diboson produc- 
tion in the same final states used for the SM Higgs 
searches are carried out. The high mass analyses measure 
pp — > W cross sections, while the low mass analyses tar- 
get VZ(— > 66) production. The data sample, reconstruc- 
tion, process modeling, uncertainties, and sub-channel 
divisions are identical to those of the SM Higgs boson 
searches. However, discriminant functions are trained to 
distinguish the contributions of SM diboson production 
from those of other backgrounds, and potential contri- 
butions from Higgs boson production are not considered. 
By way of illustration, below, we focus on VZ production. 

The NLO SM cross section for VZ production times 
the branching fraction of Z -> 66 is 0.68 ± 0.05 pb 0, 
[72j| . This is about six times larger than the 0.12 ± 
0.01 pb [2(| [26| cross section times branching fraction 
of H{-+ 66) U for a 125 GeV/c 2 SM Higgs boson, but the 
associated background is larger, due to the distribution 
of the dijet invariant mass in the U+jets events. WW 
production is considered as background. The measured 
cross section, using the MVA discriminants, for VZ is 
3.0 ± 0.6 (stat) ±0.7 fsyst) pb whereas the SM predic- 
tion is 4.4 ± 0.3 pb [34[. The combined background- 
subtracted dijet-mass distribution for the VZ analysis 
is shown in Fig. [5] for illustration. The VZ signal and 
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FIG. 2: (color online). Background-subtracted distribution 
of the reconstructed dijet mass, summed over CDF and DO's 
channels contributing to the VZ analysis. The VZ signal and 
the background contributions are fit to the data, and the fit- 
ted background is subtracted. The fitted VZ and expected 
SM Higgs (niH = 125 GeV/c 2 ) contributions are shown with 
filled histograms. The error bars shown on the data points 
correspond in each bin to the square root of the sum of the 
expected signal and background yields. 

the background contributions are fit to the data, and 
the fitted background is then subtracted. Also shown is 
the contribution expected from a SM Higgs boson with 
m,H = 125 GeV/c 2 . The W boson cross sections mea- 
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FIG. 3: (color online). Background-subtracted distribu- 
tion of the discriminant histograms, summed for bins with 
similar signal-to-background ratio (s/b) over all contribut- 
ing Higgs boson search channels from CDF and DO, for 
itlh = 125 GeV/c 2 . The background is fit to the data, and 
the uncertainty on the background, shown with the unfilled 
histogram, is after the fit. The signal model, scaled to the SM 
expectation, is shown with a filled histogram. The error bars 
shown on the data points correspond in each bin to the square 
root of the sum of expected signal and background yields. 
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FIG. 4: (color online). The log- likelihood ratio LLR as a 
function of Higgs boson mass for all of CDF and DO's SM 
Higgs boson searches in all decay modes combined. The solid 
line shows the observed LLR values, the dark long-dashed 
line shows the median expectation assuming no Higgs boson 
signal is present, and the dark- and light-shaded bands cor- 
respond, respectively, to the regions encompassing one and 
two s.d. fluctuations around the background-only expecta- 
tion. The red long-dashed line shows the median expectation 
assuming a SM Higgs boson signal is present at each value 
of rns in turn. The blue short-dashed line shows the median 
expected LLR assuming the SM Higgs boson is present at 
m H = 125 GeV/c 2 . 



sured by the high mass analyses are likewise in good 
agreement with SM predictions [H, Ell • 



B. Higgs boson combination using all decay modes 

For the search for the Higgs boson, the results pro- 
duced by the multivariate analyses can be visualized 
by combining the histograms of the final discrimi- 
nants, adding the contents of bins with similar signal- 
to-background ratio (s/b) as shown in Fig. Q] Figured] 
shows the signal expectation and the data with the back- 
ground subtracted, as a function of the s/b of the col- 
lected bins, for the combined search for a Higgs boson 
with mass mjj = 125 GeV/c 2 . The background model 
is fit to the data, allowing the nuisance parameters to 
vary within their constraints. The uncertainties on the 
background predictions in each bin are those after the fit. 
An excess of events in the highest s/b bins relative to the 
background-only expectation is observed. 

Figure @] displays the LLR distributions for the com- 
bined analyses as functions of m# . Included are the me- 
dian of the LLR distributions for the background-only hy- 
pothesis (LLRb), the signal-plus-background hypothesis 
(LLR s+ b), and the observed value for the data (LLR b s ). 
For mass hypotheses of 95 GeV/c 2 and less, fewer chan- 
nels are available for combination, giving rise to the be- 
havior of the limits shown. The shaded bands represent 
the one and two s.d. departures for LLR;, centered on the 



median. These results are listed in Table IIIII The sep- 
aration between the medians of the LLRb and LLR s+ b 
distributions provides a measure of the discriminating 
power of the search. The widths of the one- and two- 
s.d. LLRf, bands indicate the width of the LLRb distri- 
bution, assuming no signal and that fluctuations origi- 
nate from statistical fluctuations and systematic effects 
only. The value of LLR Q b s relative to LLR s+ b and LLRb 
indicates whether the data distribution more closely re- 
sembles the distributions expected if a signal is present 
(i.e., the LLR s+ b distribution, which is negative by con- 
struction) or only background is present. The signifi- 
cance of departures of LLR Q b s from LLRb can be evalu- 
ated by the width of the LLRb bands. The separation 
of the median signal-plus-background and background- 
only hypotheses is about two s.d., or greater, for Higgs 
boson masses up to ps 185 GeV/c 2 . The data are con- 
sistent with the background-only hypothesis (the black 
dashed line) at masses smaller than pa 110 GeV/c 2 and 
above approximately 145 GeV/c 2 . A slight excess is 
seen above approximately 195 GeV/c 2 , where our abil- 
ity to separate the two hypotheses is limited. For mjj 
from 115 to 140 GeV/c 2 , an excess above two s.d. in 
the data with respect to the SM background expectation 
has an amplitude consistent with the expectation for a 
standard model Higgs boson (dashed red line). Addi- 
tionally, the LLR curve under the hypothesis that a SM 
Higgs boson is present with m# = 125 GeV/c 2 is shown. 
This signal-injectcd-LLR curve has a similar shape to 



13 



the observed one. While the search for a 125 GeV/c 2 
Higgs boson is optimized to find a Higgs boson of that 
mass, the excess of events over the SM background es- 
timates also affects the results of Higgs boson searches 
at other masses. Nearby masses are the most affected, 
but the expected presence of H — > W + W~ decays for a 
125 GcV/c 2 Higgs boson implies a small expected excess 
in the H — > W + W~ searches at all masses due to the 
poor reconstructed mass resolution in this final state. 
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FIG. 6: (color online). The best-fit signal cross section ex- 
pressed as a ratio to the SM cross section as a function of 
Higgs boson mass for all of CDF and DO's SM Higgs boson 
searches in all decay modes combined. The dark- and light- 
shaded bands show the one and two s.d. uncertainty ranges on 
the fitted signal, respectively. Also shown with blue lines are 
the median fitted cross sections expected for a SM Higgs bo- 
son with uih = 125 GeV/c 2 at signal strengths of 1.0 times 
(short-dashed) and 1.5 times (long-dashed) the SM predic- 
tion. 



FIG. 5: (color online). Observed and median expected (for 
the background-only hypothesis) 95% C.L. Bayesian upper 
production limits expressed as multiples of the SM cross sec- 
tion as a function of Higgs boson mass for the combined CDF 
and DO searches in all decay modes. The dark- and light- 
shaded bands indicate, respectively, the one and two s.d prob- 
ability regions in which the limits are expected to fluctuate in 
the absence of signal. The blue short-dashed line shows me- 
dian expected limits assuming the SM Higgs boson is present 
at m H = 125 GeV/c 2 . 

The upper limit on SM Higgs boson production as a 
function of tuh is extracted in the range 90-200 GeV/c 2 
in terms of Rg]? s , the ratio of the observed limit to the pre- 
dicted SM rate. The ratios of the 95% C.L. expected and 
observed limit to the SM cross section using the Bayesian 
method are shown in Fig. [5] for the combined CDF and 
DO analyses. The observed and median-expected ratios 
are listed for the tested Higgs boson masses in Table |V] 
as obtained by the Bayesian and the CL S methods. 

Intersections of piecewise linear interpolations of the 
observed and expected rate limits with the SM=1 line are 
used to quote ranges of Higgs boson masses that are ex- 
cluded and that are expected to be excluded. The regions 
of Higgs boson masses excluded at the 95% C.L. are 90 
<m H < 109 GeV/c 2 and 149 < m H < 182 GcV/c 2 . The 
expected exclusion regions are 90 < mn < 120 GeV/c 2 
and 140 <m H < 184 GeV/c 2 . 

The observed excess for rnn from 115 to 140 GeV/c 2 
is driven by an excess of data events with respect to the 
background predictions in the most sensitive bins of the 
discriminant distributions, favoring the hypothesis that 



a signal is present. To characterize the compatibility of 
this excess with the signal-plus-background hypothesis, 
the best-fit rate cross section, i? fit , is computed, and 
shown in Fig. [6] The measured signal strength is within 
1 s.d. of the expectation for a SM Higgs boson in the 
range 115 < tuh < 140 GeV/c 2 , with maximal strength 
between 120 GeV/c 2 and 125 GcV/c 2 . At 125 GeV/c 2 , 
R Rt = 1-44+8& 

The significance of the excess in the data over the 
background prediction is computed at each hypothesized 
Higgs boson mass by calculating the local p-value un- 
der the background-only hypothesis using R m as the test 
statistic. This p-valuc expresses the probability to ob- 
tain the value of R observed in the data or larger, as- 
suming a signal is absent. These p- values are shown in 
Fig. [7] along with the expected p- values assuming a SM 
signal is present, separately for each value of m#. The 
median expected p-values assuming the SM Higgs bo- 
son is present with m#=125 GeV/c 2 for signal strengths 
of 1.0 and 1.5 times the SM prediction are also shown. 
The maximum local significance is at mn = 125 GeV/c 2 
and corresponds to 3.1 standard deviations. The fluctu- 
ations seen in the observed p-value as a function of the 
tested run result from excesses seen in different search 
channels, as well as from point-to-point fluctuations due 
to the separate discriminants at each ran, and are dis- 
cussed in more detail below. The width of the dip in the 
observed p-values from 115 to 140 GcV/c 2 is consistent 
with the resolution of the combination of the H — > bb and 
H —t W + W~ channels, as illustrated by the injected 
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FIG. 7: (color online). The solid black line shows the back- 
ground p-value as a function of ms for all of CDF and DO's 
SM Higgs boson searches in all decay modes combined. The 
dotted black line shows the median expected values assum- 
ing a SM signal is present, evaluated separately at each tuh- 
The associated dark- and light-shaded bands indicate the 
one and two s.d. fluctuations of possible experimental out- 
comes under this scenario. The blue lines show the median 
expected p-values assuming the SM Higgs boson is present 
with m_ff = 125 GeV/c 2 at signal strengths of 1.0 times (short- 
dashed) and 1.5 times (long-dashed) the SM prediction. 



FIG. 8: (color online). The solid black line shows the signal- 
plus-background p-value as a function of mH for all of CDF 
and DO's SM Higgs boson searches in all decay modes com- 
bined. The dotted black line shows the median expected val- 
ues assuming no SM signal is present, evaluated separately 
at each uih- The associated dark and light-shaded bands 
indicate the one and two s.d. fluctuations of possible exper- 
imental outcomes under this scenario. The blue lines show 
the median expected p-values assuming the SM Higgs bo- 
son is present with m_ff=125 GeV/c 2 at signal strengths of 
1.0 times (short-dashed) and 1.5 times (long-dashed) the SM 
prediction. 



signal curves in Fig. [7] The effective resolution of this 
search comes from two independent sources of informa- 
tion. The reconstructed candidate masses help constrain 
to//, but more importantly, the expected cross sections 
times the relevant branching ratios for the H — > bb and 
H — > W + W~ channels are strong functions of to,// in the 
SM. The observed excess in the H -> bb channels coupled 
with the slight excess in the H — > W + W~ channels de- 
termine the shape of the observed p-value as a function 
of to//. 

Figure [8] shows the quantity CL s+ b, corresponding to 
the p-value for the signal-plus-background hypothesis. 
The observed value, along with the expected p-values as- 
suming a signal is absent, are shown separately for each 
value of to//. The median expected p-values assuming the 
SM Higgs boson is present with to//=125 GeV/c 2 for sig- 
nal strengths of 1.0 and 1.5 times the SM prediction are 
also shown. In the mass region from 115 to 140 GeV/c 2 
the observed values above 50% indicate a high level of 
consistency with the signal-plus-background hypothesis. 

We also separate CDF and DO's searches into combina- 
tions focusing on the H bb, H — >■ W + W~ , H — > 77, 
and H —> t + t~ decay modes, and these are discussed 
in the following sections. 

C. H ^ bb Decay Mode 

Below 130 GeV/c 2 , the H — > bb searches contribute 
the majority of the search sensitivity. The WH — > Ivbb, 
ZH —> vvbb, and ZH l + t~bb channels from both ex- 



periments are included in this sub-combination. Two of 
the six contributing channels were updated for this sub- 
combination compared with that reported in Ref. [l2j]. 
The CDF ZH — > vvbb [73| analysis was updated to use a 
more powerful MVA b-tagging algorithm [44| along with 
changes to the kinematic selections. The assignment of 
correlated systematic uncertainties between channels was 
updated in the DO WH -> Ivbb analysis [II]. The ob- 
served LLR distribution is shown in Fig. HI along with its 
expected values under the background-only and signal- 
plus-background hypotheses. The hypotheses that a SM 
Higgs boson is present with to// = 125 GeV/c 2 for signal 
strengths of 1.0 and 1.5 times the SM prediction are also 
given. The LLR values as a function of Higgs boson mass 
are listed in Table ED 

We multiply the best-fit rate cross section, R , for this 
sub-combination by the SM prediction for the associated- 
production cross section times the decay branching ratio 
ipwn + &zh) x B{H — > bb), to obtain the observed value 
for this quantity. We show the fitted (own + &zh) x 
B(H — > bb) as a function of to//, along with the SM 
prediction, in Fig. [TO] The figure also shows the expected 
cross section fits for each to//, assuming that the SM 
Higgs boson with to// = 125 GeV/c 2 is present, both at 
the rate predicted by the SM, and also at a multiple of 
1.5 times that of the SM. The best-fit rate corresponds to 
(own + <?zh) x B(H -> bb) = 0.19±g;g| (stat + syst) pb. 
The shift in this result compared with the value of 0.23 ± 
0.09 (stat + syst) pb obtained previously |12l| is due to the 
updated ZH —> vvbb analysis from CDF [50l . [73| , and 
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FIG. 9: (color online). The log-likelihood ratio LLR as a 
function of Higgs boson mass from the combination of CDF 
and DO's SM Higgs boson searches focusing on the H — > 66 de- 
cay mode. The solid line shows the observed LLR values, the 
dark short-dashed line shows the median expectation assum- 
ing no Higgs boson signal is present, and the dark- and light- 
shaded bands correspond, respectively, to the regions encom- 
passing one and two s.d. fluctuations around the background- 
only expectation. The red dash-dot line shows the median 
expectation assuming a SM Higgs boson signal is present at 
each value of ttih in turn. The blue lines show the median 
expected LLR assuming the SM Higgs boson is present at 
itlh = 125 GeV/c 2 with signal strengths of 1.0 times (short- 
dashed) and 1.5 times (long-dashed) the SM prediction. 
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FIG. 10: (color online). The best-fit signal cross section 
times branching ratio (pwH + cfzh) X B(H — >■ 66) as a func- 
tion of Higgs boson mass from the combination of CDF and 
DO's SM Higgs boson searches focusing on the H — > bb decay 
mode. The dark- and light-shaded bands show the one and 
two s.d. uncertainty ranges on the fitted signal, respectively. 
Also shown with blue lines are the median fitted cross sec- 
tions expected for a SM Higgs boson with rnn = 125 GeV/c 2 
at signal strengths of 1.0 times (short-dashed) and 1.5 times 
(long-dashed) the SM prediction. The SM prediction is shown 
as the smooth, falling curve where the narrow band indicates 
the theoretical uncertainty. 



corresponds to a change in the central value of 0.5 times 
the total uncertainty. For ran = 125 GeV/c 2 , the SM 
predicts (a WH + a ZH ) x B(H -> bb) = 0.12± 0.01 pb. 



D. H -y W+W Decay Mode 

Above 130 GeV/c 2 , the H W+W- channels con- 
tribute the majority of the search sensitivity. We com- 
bine all H — > W + W~ searches from CDF and DO, incor- 
porating potential signal contributions from gluon-gluon 
fusion, WH, ZH, and vector boson fusion production. 
Approximately 75% of the signal comes from the gluon- 
gluon fusion process, 20% from associated production, 
and 5% from the VBF process. The LLR distributions 
are shown in Fig. [TT] and the values as a function of Higgs 
boson mass are listed in Table IIVI The data present a 
one to two s.d. excess in the region from 115 to 140 
GeV/c 2 where there is some separation between the two 
hypotheses. An excess is also seen in the searches for 
Higgs bosons with mass mjj > 195 GeV/c 2 , as men- 
tioned in Section IVIII B[ but the sensitivity to the SM 
Higgs boson is not as large at these masses as it is at 
lower masses. Figure [T2l shows the best-fit cross section 
for the combined H — > W+W~ searches, normalized to 
the SM prediction, as a function of tuh , along with the 
expectations assuming the SM Higgs boson is present at 
m,H = 125 GeV/c 2 for signal strengths of 1.0 and 1.5 
times the SM prediction. 



E. H — > 77 Decay Mode 

We also separately combine CDF and DO's searches 
focusing on the H — > 77 decay mode and display the re- 
sulting upper limits on the production cross section times 
the decay branching ratio normalized to the SM predic- 
tion in Fig. 1131 An excess of approximately two s.d. is 
seen in these searches at mjj = 125 GeV/c 2 , but its con- 
tributions to the fully combined SM cross section and 
limit are small due to the low expected signal yield in this 
channel. However, the observed excess in the H — > 77 
search channel has a visible impact on Higgs boson cou- 
pling constraints as described in Section IVlII Gl 

F. H — >■ tt Decay Mode 

We also separately combine CDF and DO's searches fo- 
cusing on the H — > t + t~ decay mode and display the re- 
sulting upper limits on the production cross section times 
the decay branching ratio normalized to the SM predic- 
tion in Fig. [lH 
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FIG. 11: (color online). The log-likelihood ratio LLR as 
a function of Higgs boson mass from the combination of 
CDF and DO's SM Higgs boson searches focusing on the 
H — > W + W~ decay mode. The solid line shows the observed 
LLR values, the dark short-dashed line shows the median ex- 
pectation assuming no Higgs boson signal is present, and the 
dark- and light-shaded bands correspond, respectively, to the 
regions encompassing one and two s.d. fluctuations around the 
background-only expectation. The red dash-dot line shows 
the median expectation assuming a SM Higgs boson signal 
is present at each value of tuh in turn. The blue lines show 
the median expected LLR assuming the SM Higgs boson is 
present at inn = 125 GeV/c 2 with signal strengths of 1.0 
times (short-dashed) and 1.5 times (long-dashed) the SM pre- 
diction. 
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FIG. 13: (color online). Observed and median expected (for 
the background-only hypothesis) 95% C.L. Bayesian upper 
production limits expressed as multiples of the SM cross sec- 
tion as a function of Higgs boson mass from the combination 
of CDF and DO's SM Higgs boson searches focusing on the 
H — > 77 decay mode. The dark- and light-shaded bands in- 
dicate, respectively, the one and two s.d. probability regions 
in which the limits are expected to fluctuate in the absence 
of signal. 
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FIG. 12: (color online). The best-fit signal cross section 
expressed as a ratio to the SM cross section as a function 
of Higgs boson mass from the combination of CDF and DO's 
SM Higgs boson searches focusing on the H — > W + W~ decay 
mode. The dark- and light-shaded bands show the one and 
two s.d. uncertainty ranges on the fitted signal, respectively. 
Also shown with blue lines are the median fitted cross sec- 
tions expected for a SM Higgs boson with rnn = 125 GeV/c 2 
at signal strengths of 1.0 times (short-dashed) and 1.5 times 
(long-dashed) the SM prediction. 
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FIG. 14: (color online). Observed and median expected (for 
the background-only hypothesis) 95% C.L. Bayesian upper 
production limits expressed as multiples of the SM cross sec- 
tion as a function of Higgs boson mass from the combination 
of CDF and DO's SM Higgs boson searches focusing on the 
H — > t + t~ decay mode. The dark- and light-shaded bands 
indicate, respectively, the one and two s.d. probability regions 
in which the limits are expected to fluctuate in the absence 
of signal. 
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TABLE III: Log-likelihood ratio (LLR) values obtained from the combination of all of CDF and DO's Higgs boson search 
channels using the CL S method. 



m H (GeV/c 2 ) 


LLR b s 


LLR s+ 5 


LLR- 2ff 


LLR^ lCT 


LLR;, 


LLR+ lCT 


llr+ 2ct 


90 


17.02 


-7.24 


17.31 


12.08 


6.84 


1.61 


-3.62 


95 


13.07 


-5.96 


15.21 


10.44 


5.68 


0.91 


-3.85 


100 


8.39 


-7.44 


17.73 


12.40 


7.08 


1.76 


-3.56 


105 


3.62 


-6.69 


16.38 


11.35 


6.32 


1.29 


-3.74 


110 


2.53 


-5.73 


14.79 


10.12 


5.45 


0.78 


-3.89 


115 


-3.67 


-4.81 


13.17 


8.88 


4.59 


0.31 


-3.98 


120 


-8.44 


-4.09 


11.76 


7.82 


3.88 


-0.06 


-4.00 


125 


-7.72 


-3.52 


10.76 


7.07 


3.39 


-0.29 


-3.97 


130 


-3.74 


-3.30 


10.31 


6.74 


3.18 


-0.39 


-3.95 


135 


-4.81 


-3.64 


10.89 


7.17 


3.45 


-0.26 


-3.98 


140 


-5.08 


-4.09 


11.72 


7.79 


3.86 


-0.07 


-4.00 


145 


0.20 


-5.07 


13.35 


9.02 


4.69 


0.36 


-3.97 


150 


3.72 


-6.68 


15.87 


10.95 


6.04 


1.12 


-3.79 


155 


8.44 


-8.80 


18.72 


13.18 


7.65 


2.12 


-3.41 


160 


13.45 


-15.25 


26.04 


19.08 


12.12 


5.15 


-1.81 




1 7 


— 1 7 81 
n .01 


98 7(\ 


91 ^l 


1 3 87 

lO.OI 


(\ 49 


—i rn 


170 


10.93 


-12.26 


22.87 


16.50 


10.13 


3.77 


-2.60 


175 


7.33 


-8.77 


18.50 


13.02 


7.53 


2.04 


-3.45 


180 


4.86 


-6.17 


14.87 


10.18 


5.50 


0.81 


-3.88 


185 


2.14 


-3.92 


11.23 


7.42 


3.62 


-0.19 


-3.99 


190 


-0.99 


-2.61 


8.73 


5.60 


2.46 


-0.68 


-3.81 


195 


-2.83 


-1.98 


7.34 


4.60 


1.87 


-0.87 


-3.60 


200 


-2.50 


-1.53 


6.29 


3.88 


1.46 


-0.96 


-3.37 


TABLE IV: Log 


;-likelihood ratio (LLR) values obtained from the combination i 


Df CDF and DO's 


Higgs boson search channels 


focusing on the H — > W + W decay mode using 


; the CL S method. 










m H (GcV/c 2 ) 


LLR bs 


LLRs+t, 


LLR" 2a 


LLR~ la 


LLR& 


LLR+ lff 


LLR+ 2ct 


100 


-0.42 


-0.10 


1.36 


0.73 


0.10 


-0.53 


-1.16 


105 


-0.72 


-0.10 


1.34 


0.72 


0.10 


-0.53 


-1.15 


110 


-1.07 


-0.14 


1.62 


0.88 


0.14 


-0.60 


-1.34 


115 


-0.80 


-0.24 


2.18 


1.21 


0.24 


-0.74 


-1.71 


120 


-0.98 


-0.49 


3.26 


1.87 


0.48 


-0.91 


-2.30 


125 


-1.69 


-0.94 


4.73 


2.82 


0.91 


-1.00 


-2.91 


130 


-0.59 


-1.57 


6.44 


3.98 


1.52 


-0.95 


-3.41 


135 


-1.11 


-2.40 


8.36 


5.33 


2.30 


-0.73 


-3.77 


140 


-3.38 


-3.24 


10.09 


6.58 


3.08 


-0.43 


-3.94 


145 


1.19 


-4.42 


12.23 


8.17 


4.12 


0.06 


-4.00 


150 


3.43 


-6.13 


15.01 


10.29 


5.57 


0.85 


-3.87 


155 


8.05 


-8.59 


18.45 


12.97 


7.50 


2.02 


-3.46 


160 


13.27 


-15.15 


25.92 


18.98 


12.04 


5.10 


-1.84 


165 


17.55 


-17.75 


28.69 


21.25 


13.82 


6.38 


-1.05 


170 


11.19 


-12.21 


22.80 


16.45 


10.09 


3.74 


-2.61 


175 


7.28 


-8.72 


18.44 


12.96 


7.49 


2.02 


-3.46 


180 


4.63 


-6.12 


14.80 


10.13 


5.46 


0.78 


-3.89 


185 


1.56 


-3.83 


11.05 


7.29 


3.53 


-0.23 


-3.99 


190 


-1.39 


-2.50 


8.51 


5.44 


2.36 


-0.71 


-3.79 


195 


-3.24 


-1.88 


7.12 


4.45 


1.78 


-0.89 


-3.56 


200 


-3.23 


-1.45 


6.08 


3.73 


1.38 


-0.97 


-3.32 
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TABLE V: Ratios of observed and median expected (for the 
background-only hypothesis) 95% C.L. upper production lim- 
its to the SM cross section as a function of the Higgs boson 
mass for the combined CDF and DO searches in all decay 
modes, obtained using the Bayesian and CL S methods. 



Baycsian 



CL S 



mu (GeV/c 2 ) 


nobs 
^95 


1L 95 


Dobs 
^-95 


-"-95 


90 


0.37 


0.74 


0.39 


0.74 


95 


0.48 


0.80 


0.49 


0.81 


100 


0.62 


0.72 


0.62 


0.73 


105 


0.89 


0.77 


0.93 


0.77 


110 


1.02 


0.82 


1.03 


0.83 


115 


1.63 


0.90 


1.67 


0.91 


120 


2.33 


1.00 


2.40 


0.99 


125 


2.44 


1.06 


2.62 


1.07 


130 


2.13 


1.11 


2.10 


1.10 


135 


2.03 


1.04 


2.12 


1.06 


140 


2.10 


1.01 


2.08 


1.00 






n 88 

u.oo 


1 9Q 


n on 


150 


0.94 


0.79 


0.91 


0.78 


155 


0.64 


0.69 


0.62 


0.68 


160 


0.46 


0.51 


0.45 


0.51 


165 


0.37 


0.47 


0.36 


0.47 


170 


0.54 


0.57 


0.53 


0.57 


175 


0.71 


0.68 


0.68 


0.68 


180 


0.87 


0.81 


0.86 


0.82 


185 


1.20 


1.02 


1.18 


1.04 


190 


1.86 


1.29 


1.86 


1.27 


195 


2.74 


1.44 


2.64 


1.48 


200 


3.07 


1.66 


2.97 


1.67 



G. Compatibility of the Excess with the SM Higgs 
Boson Hypothesis 

The best-fit rate parameters, R Rt , for the full com- 
bination of all channels and the combinations of chan- 
nels focusing on the H — > W + W~,H — > bb, H — > 77, 
and H — > t + t~ decay modes [74j] are listed in Ta- 
ble IVIII as a function of Higgs boson mass over the range 
115 < m,H < 140 GeV/c 2 , where the combined result 
has sensitivity to a signal and a clear excess exists. For 



m H = 125 GeV/c 2 , we obtain R ttt = 1.44±g;|| using all 
decay modes. 

Figure [15] shows the contribution of the four combina- 
tions for the different decay modes to the best-fit signal 
cross section for inn = 125 GeV/c 2 . The results arc 
consistent with each other, with the full combination, 
and with the production of the SM Higgs boson at that 
mass. Figure \W\ shows the posterior probability densi- 
ties obtained for the cross section scale factors from the 
H ->• bb, H -> W + W~, H -> 77, and H -> t+t~ com- 
binations. 

The Higgs boson is expected to couple more strongly 
to more massive particles than to less massive ones, and 
thus may provide sensitivity to non-SM particles whose 
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FIG. 15: (color online) Best-fit values of R = {a x S)/SM 
using the Bayesian method for the combinations of CDF 
and DO's Higgs boson search channels focusing on the H — > 
W + W~ , H — > bb, H 77, and H — > t + t~ decay modes for 
a Higgs boson mass of 125 GeV/c 2 . The shaded band corre- 
sponds to the one s.d. uncertainty on the best-fit value of -R 
for all SM Higgs boson decay modes combined. 



interactions become more relevant at higher energies. It 
is important therefore to study in detail the properties 
of the new particle. The channel- by- channel values of 
R = (a x B)/SM provide useful constraints on the possi- 
ble couplings of the particle [75j , but their interpretation 
is ambiguous because signal contributions from multiple 
sources are simultaneously accepted by each sub-channel. 
For example, the ZH — ► vvbb channels have sensitivity 
to both the WH and ZH production modes, and the 
H W + W~ searches are sensitive to gluon-gluon fu- 
sion, WH, ZH, and VBF in different mixtures within 
independent sub-channels characterized by the number 
of reconstructed jets. 

Most of the searches conducted at the Tevatron are 
sensitive to the product of fermion and boson coupling 
strengths. In the VH — > Vbb searches, the production 
depends on the coupling of the Higgs boson to the weak 
vector bosons, while the decay is to fermions. In the 
gg —> H — >• W + W~ searches, the production is domi- 
nated by the Higgs boson couplings to fermions via the 
quark loop processes, but the decay is to bosons. A large 
enhancement of the Higgs boson's couplings to fermions 
can thus be masked by a small coupling to bosons, and 
vice versa, as shown in Fig. 2 of Ref. [75[. However, 
other less-sensitive channels included in this combina- 
tion provide additional constraints. The same-sign di- 
lcpton searches, the tri-lepton searches, and some of the 
searches with tau leptons as decay products of W bosons 
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TABLE VI: Log-likelihood ratio (LLR) values obtained from the combination of CDF and DO's Higgs boson search channels 
focusing on the H — >• bb decay mode using the CL S method. 



rriH (CjcV / c ) 


i i T) 

LLKobs 




T T D —2(7 

LL\x h 


I I T> — 1(7 

LL\x b 


I 1 I > 

LLiXb 


LLR b 


LL\x h 


90 


1 +7 AO 

17.02 


+7 O A 
— 7.24: 


17.31 


12.08 


f O A 

6.84 


1.61 


—3.62 


95 


1 O AT 

13.07 


— 5.96 


15.21 


10.44 


5.68 


a a i 
0.91 


— 3.85 


100 


9.50 


—6.71 


16.41 


11.38 


6.34 


1.30 


—3.73 


105 


6.09 


—6.00 


15.12 


1 A OO 

10.38 


5.63 


A O O 

0.88 


O Of 

—3.86 


110 


2.21 


-5.06 


13.52 


9.15 


4.78 


0.41 


-3.97 


115 


-2.84 


-4.12 


11.83 


7.87 


3.91 


-0.04 


-4.00 


120 


-6.05 


-3.16 


9.91 


6.45 


2.99 


-0.47 


-3.93 


125 


-5.05 


-2.16 


7.88 


4.99 


2.09 


-0.80 


-3.69 


130 


-3.53 


-1.34 


5.89 


3.60 


1.31 


-0.98 


-3.27 


135 


-4.60 


-0.80 


4.28 


2.52 


0.77 


-0.98 


-2.74 


140 


-1.77 


-0.41 


2.93 


1.67 


0.40 


-0.87 


-2.13 


145 


-0.86 


-0.20 


1.95 


1.07 


0.19 


-0.68 


-1.56 


150 


-0.31 


-0.08 


1.23 


0.66 


0.08 


-0.49 


-1.07 



TABLE VII: Best-fit values of R = (a x S)/SM using the Bayesian method for all SM Higgs boson decay modes combined 
and the combinations of CDF and DO's Higgs boson search channels focusing on the H — > W + W~ , H — » bb, H — » 77, and 
H — > t + t~ decay modes as a function of Higgs boson mass over the range 115 < ran < 140 GeV/c 2 . The quoted uncertainties 
bound the smallest interval containing 68% of the integral of the posterior probability density. 



m H (GeV/c 2 ) 


115 


120 


125 


130 


135 


140 


ifet(SM) 


0.8218:8 


1 49+0.53 


1 44+0-59 


-, 10+0. 60 
i ' lo -0.60 


n qq+0-58 


1 ic+0.57 
1 - io -0.52 


R R t{H -> W+W-) 


n 00+1-65 
z - zz -1.59 


1 t:q+1.20 


Q4+0-85 
u - y ^-0.83 


4Q+0-69 
u -^ y -0.63 


n M+0-53 
u - O4 -0.52 


0.97l°;| 


R m (H^bb) 


a 7 9+0.47 


1 26+ 062 

1 - ZO -0.55 


1.59t° f 2 


1 (59+0.91 


2 ^9 + 1.22 
Z - DZ -1.21 


9q+l-61 
° — 1.74 


Rat(H -> 77) 


RET+ 2 - 66 

U - DO -0.54 


c 04+3.20 


5.97H:» 


3.17t 2 ; 6 8 ? 


u - uu -o.oo 


01 +3.30 
J - 01 -3.13 


R m {H^T+T-) 


1 70+ 2 - 20 


no +2 - 22 


i-68i?:l 


O-OOt 2 ;^ 


oo+ 2 - 83 
u - uu -o.oo 


-1 9C+2.62 
1 - zo -1.15 



are primarily sensitive to VH — > VW + W~ , an entirely 
bosonic process, although their results are customarily 
reported in combination with the other H — > W + W~ 
searches. The searches for tiH — > ttbb provide constraints 
on the fermion couplings with minimal masking from the 
bosonic couplings. 

We follow the notation of Ref. [76| and introduce mul- 
tiplicative scaling factors for the coupling of the Higgs 
boson to fermions (k/) and either to W bosons (kw) 
and Z bosons (kz) or more generically to vector bosons 
(kv)- We then search for deviations from the expected 
SM values of /tj = 1. 

The first test assumes ran = 125 GeV/c 2 , based 
on the ATLAS and CMS observations, and fits for the 
H —> W + W~ coupling, holding all other couplings fixed 
to their SM values. This test corresponds to holding 
the values of k,z = ft/ = 1, while varying nyy- At each 
value of Kw, we recompute the predicted cross sections 
and decay branching ratios, as described in Ref. [76|. We 
assume a uniform prior density in kw, and show the pos- 



terior probability distribution in Fig. 1171 A negative sign 
of kw is preferred by the Tevatron data due to the excess 
seen in the H — > 77 searches. In the SM, this process 
proceeds at lowest order via a VF-boson loop or a quark 
loop (dominated by the top quark), with destructive in- 
terference between the two contributions [77|, as given 
by T(H -> 77) = T(H -> 77)5^/ X 1 1.28 k v - 0.21 n f \ 2 . 
If the sign of the H — > W + W~ coupling is negative, 
then this interference becomes constructive, allowing for 
a larger prediction of the yield. We obtain a best-fit value 
of Ky/ = —1.27. Our procedure for finding the smallest 
set of intervals that contain 68% of the integral of the 
posterior results in two intervals, —1.56 < % < —0.81 
and 1.04 < Kw < 1-51. We perform a similar test for 
kz, assuming kw = kj = 1. The resulting posterior 
probability density is shown in Fig. [T5] The Higgs boson 
searches at the Tevatron are sensitive almost exclusively 
to the square of Kz, and thus the posterior density is 
nearly symmetric in positive and negative couplings. The 
best-fit values are k z = 1.05±g;|| and k z = -lMlio'.ll- 
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FIG. 16: (color online), (a) Posterior probability densities 
for R = (a x S)/SM using the Bayesian method from the 
combinations of CDF and DO's Higgs boson search channels 
focusing on the H — ¥ W + W~, H — > bb, H — > 77, and H — > 
t + t~ decay modes and for all SM Higgs boson decay modes 
combined. The same curves are shown on an expanded scale 
in (b). 

Finally, we perform a similar test for Kf, the common 
scale factor on the Higgs boson couplings to fermions, 
holding kw = Kz = 1- The resulting posterior prob- 
ability density is shown in Fig. 1191 An asymmetry is 
seen in this distribution, due again to the outcome in 
the H — > 77 channels. We obtain a best-fit value of 
Kf = — 2.64+^3Q. The large magnitude of the fitted value 
is due to the excesses seen in the H — > bb and H — > 77 
searches. 

We then allow both kw and kz to vary independently, 
also allowing to vary by integrating the likelihood 
function times a uniform prior in kj over negative and 
positive values. The resulting areas in the (kw, Kz) plane 
preferred by the Tevatron data are shown in Fig. [2111 
While we allow either coupling scale factor to be negative, 
only two quadrants are shown in Fig. [20] due to an overall 
sign ambiguity. The point (kw, k z) = (0,0) corresponds 
to no Higgs boson production or decay in the most sen- 
sitive search modes at the Tevatron and is excluded at 
more than the 95% C.L. due to the Higgs-boson-like sig- 



FIG. 17: (color online). Posterior probability density for 
Kw from the combination of Tevatron searches for a SM-like 
Higgs boson with mn = 125 GeV/c 2 . The couplings of the 
Higgs boson to fermions and to the Z boson are assumed to 
be as predicted by the SM. The values that maximize local 
posterior probability densities are shown with dashed lines, 
and the 68% and 95% C.L. intervals are indicated with the 
dark- and light-shaded regions, respectively. The predicted 
SM value of kw is indicated by the solid line. 
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FIG. 18: (color online). Posterior probability density for 
Kz from the combination of Tevatron searches for a SM-like 
Higgs boson with ran = 125 GeV/c 2 . The couplings of the 
Higgs boson to fermions and to the W ± boson are assumed 
to be as predicted by the SM. The values that maximize local 
posterior probability densities are shown with dashed lines, 
and the 68% and 95% C.L. intervals are indicated with the 
dark- and light-shaded regions, respectively. The predicted 
SM value of kz is indicated by the solid line. 



nal in the H —¥ bb and H — > W + W channels. Our 
best-fit points are (kw,kz) = (1.25, ±0.90). 

We study the ratio Xwz = kw/kz using the same 
posterior probability density that is used in Fig. [20] 
We choose a projection onto a one dimensional vari- 
able that preserves the uniformity of the prior proba- 
bility density in the two-dimensional plane. This vari- 
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FIG. 19: (color online). Posterior probability density for 
Kf from the combination of Tevatron searches for a SM-like 
Higgs boson with mn = 125 GeV/c 2 . The couplings of the 
Higgs boson to the W ± and Z bosons are assumed to be as 
predicted by the SM. The values that maximize local posterior 
probability densities are shown with dashed lines, and the 68% 
and 95% C.L. intervals are indicated with the dark- and light- 
shaded regions, respectively. The predicted SM value of Kf is 
indicated by the solid line. 



FIG. 21: (color online). Posterior probability density for 
Owz = tan _1 («;z / r Kw), from the combination of Tevatron 
searches for a SM-like Higgs boson with ran = 125 GeV/c 2 
allowing Kf to float. The value that maximizes the poste- 
rior probability density is shown with a dashed vertical line, 
and the 68% and 95% C.L. intervals are indicated with the 
dark- and light-shaded regions, respectively. The predicted 
SM value of Owz is indicated by the solid vertical line. 




FIG. 20: (color online). Two-dimensional constraints in 
the (kw,kz) plane, for the combined Tevatron searches for 
a SM-like Higgs boson with mass 125 GeV/c 2 allowing Kf 
to float. The points that maximize local posterior probability 
densities are marked with dots, and the 68% and 95% C.L. in- 
tervals are indicated with the dark- and light-shaded regions, 
respectively. The SM prediction for (kw,kz) is marked with 
a triangle. 



FIG. 22: (color online). Two-dimensional constraints in the 
(kv, Kf) plane, for the combined Tevatron searches for a SM- 
like Higgs boson with mass 125 GeV/c 2 assuming Custodial 
symmetry (\wz = 1). The points that maximize local poste- 
rior probability densities are marked with dots, and the 68% 
and 95% C.L. intervals are indicated with the dark- and light- 
shaded regions, respectively. The SM prediction for (kv, Kf) 
is marked with a triangle. 
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able is the angle with respect to the Kw axis, 6\yz — 
tan -1 (kz/kw) = tan _1 (l/Aiyz). Figure [23] shows the 
one-dimensional posterior probability density in this vari- 
able. This function is symmetric for positive and negative 
Owz- We measure \6 



wz\ 



0.68j^4i, which corresponds 



to A 



wz 



1.24+1 



.34 
0.42' 



Assuming that custodial symmetry |78j holds (Xwz = 
1), we allow both ny and Kf to vary, and show in Fig. 1221 
the regions preferred at the 68% C.L. and the 95% C.L. 
in the two-dimensional plane («y,K/). The asymme- 
try induced by the excesses in the H — > 77 searches is 
visible in this projection as well. The best-fit point is 
(Ky,Kf) = (1.05,-2.40), but a secondary maximum in 
the posterior density is seen at (Ky,K/) = (1.05,2.30), 
consistent with the SM expectation, given the large un- 
certainties. The integral of the posterior density in the 
(+,+) quadrant is 26% of the total, while the remaining 
74% of the integral of the posterior density is contained 
within the (+,-) quadrant. 



IX. RESULTS - NON-STANDARD MODEL 
INTERPRETATIONS 

The mechanism of elcctroweak symmetry breaking 
may entail a richer phenomenology than expected in the 
SM. Natural extensions include the addition of a fourth 
generation of fcrmions with masses much larger than 
those of the three known generations or models with sev- 
eral Higgs bosons or models in which the Higgs boson(s) 
may have modified couplings. We interpret our Higgs 
boson search results in models with a sequential fourth 
generation of fermions (SM4) and in the fermiophobic 
Higgs model (FHM) described below. 



A. Fourth Generation Interpretation 

With the inclusion of two additional heavy fourth- 
generation quarks in the SM4 [79], the gg — > H cou- 
pling is enhanced b y a factor of roughly three relative 
to the SM coupling [80l - l82| . The partial decay width for 
H — > gg is enhanced by the same factor as the produc- 
tion cross section. However, because the H — >■ gg decay 
is mediated by a loop amplitude, the H — > W + W~ de- 
cay continues to dominate for Higgs boson masses above 
135 GeV/c 2 . Since the expected signal yield is larger in 
the SM4 model than the SM, the sensitivity of CDF and 
DO's Higgs boson searches extends to higher masses. For 
this reason, the upper end of the search range for the rel- 
evant channels is raised to 300 GeV/c 2 for interpretations 
associated with this model. 

Two scenarios for the masses of the fourth-generation 
fermions are considered. In the first, the low-mass sce- 
nario, we set the mass of the fourth-generation neutrino 
to„4 = 80 GeV/c 2 and the mass of the fourth-generation 
charged lepton = 100 GeV/c 2 , to be compatible with 



current experimental constraints [83j and have the max- 
imum impact on the Higgs boson decay branching ra- 
tios. In the second, the high-mass scenario, we set m„4 = 
iri(4 = 1 TeV/c 2 , so that the fourth-generation leptons do 
not modify the decay branching ratios of the Higgs bo- 
son relative to the SM. In both scenarios, we choose the 
masses of the quarks to be those of the second scenario 
in Ref. H {m dA = 400 GeV/c 2 and m uA = 450 GeV/c 2 ). 
The next-to-next-to-leading order (NNLO) production 
cross section calculation of Ref. (HJ is used, which is a 
modified version of the NNLO SM calculation. Previous 
interpretations of SM Higgs boson searches within the 
context of a fourth generation of fermions at the Tevatron 
excluded 131 < ran < 207 GeV/c 2 |84|. Similar searches 
have been performed by the ATLAsjH and CMS HI 
Collaborations, excluding 140 < tuh < 185 GeV/c 2 and 
144 < nriH < 207 GeV/c 2 , respectively. A more re- 
cent search by the CMS Collaboration excluded the mass 
range 110 < m H < 600 GeV/c 2 H3. 

We combine our searches for a Higgs boson in the pro- 
cesses gg — >• H — >• W + W~ and gg — >• H — >• ZZ. Limits 
on the SM4 models and on <r(gg -> H)xB{H -> W + W~) 
are derived. This result is an update of Ref. |83 |. The 
analyses are performed equivalently to the SM searches 
except that gg — >■ H production only is considered for 
the signal. The MVA classifiers are retrained accord- 
ingly and, for the specific case of the DO H —> W + W~ — > 
l+vl^v channel, the two-jet bin, which is less sensitive 
to gg — > H production, is not included. 

The branching ratios for H — > W + W~ are calculated 
using hdecay [231 modified to include fourth-generation 
fermions [8l[. To include the gg — > H — > ZZ searches, we 
assume the SM value for B(H -> W + W~)/B{H ->• ZZ). 
In setting limits on a(gg —> H ) x B(H —> W + W~), the 
gg — > H — > ZZ process is included assuming that its 
signal yield scales equivalently to that from the gg — > 
H -> W+W~ channel. 

When setting limits on a(gg — > H) xB(H — s- I / F + M /_ ), 
the theoretical uncertainty on the prediction of o~(gg — > 
H) x B(H — » IT + VF~) is not included since these lim- 
its are independent of the predictions. However, when 
setting limits on m# in the context of fourth-generation 
models, uncertainties on the theoretical predictions are 
included as described for the SM searches. 

The combined limits on a(gg H) x B(H -> W + W~) 
obtained using the Baycsian method are shown in Fig. [23] 
along with the theory predictions for fourth-generation 
models in the low- and high-mass scenarios. Limits ob- 
tained using both the Bayesian and CL S methods are 
listed as a function of Higgs boson mass in Table IVIIII 
A broad, moderate excess above the background expec- 
tation is seen for masses above 200 GeV/c 2 . 

Production limits obtained for the two SM4 scenarios 
using the Bayesian method are shown in Fig. [M] The 
limits are presented as ratios relative to SM4 low-mass 
scenario predictions as a function of the Higgs boson 
mass. In the low-mass scenario, which gives the smaller 
excluded mass range, a SM-likc Higgs boson with a mass 



23 



10 



x 

T 

s 1 



10 



Tevatron Run II, L int < 10 fb 
SM4 Higgs combination 



¥*6 



Exp. 95% C.L Limit (w/o Higgs) 
Obs. 95% C.L. Limit 
Expected ± 1 s.d. 
Expected ± 2 s.d. 
SM4 (low mass) 
SM4 (high mass) 




100 120 140 160 180 200 220 240 260 280 300 

m H (GeV/c 2 ) 



FIG. 23: (color online). Observed and median expected (for 
the background-only hypothesis) 95% C.L. Bayesian upper 
limits on the cross section times branching ratio a(gg — > 
H) x B(H -> W + W~) from the combination of CDF and 
DO's Higgs boson search channels focusing on this production 
and decay mode. The dark and light-shaded bands indicate, 
respectively, the one and two s.d. probability regions in which 
the limits are expected to fluctuate in the absence of signal. 
Theoretical predictions for SM4 in the low- and high-mass sce- 
narios are shown with blue and red lines, respectively. The 
hatched band indicates the theoretical uncertainty associated 
with the SM4 low-mass scenario. 



in the range 121-225 GeV/c 2 is excluded at the 95% C.L. 
The expected excluded mass range is 118-270 GeV/c 2 . In 
the high-mass scenario, the mass range 121-232 GeV/c 2 
is excluded, with an expected excluded mass range of 
118-290 GeV/c 2 . 



B. Fermiophobic Interpretation 

In the FHM, the lightest Higgs boson does not cou- 
ple to fermions at tree level, but aside from this one 
difference, its behavior is indistinguishable from that of 
the SM Higgs boson. In the FHM, the production of 
Higgs bosons, Hf, at hadron colliders via the process 
gg — > Hf is suppressed to a negligible rate and is ig- 
nored in the context of this interpretation. The as- 
sociated production mechanisms pp — > WHf + X and 
pp — > ZHf + X, as well as the vector-boson-fusion (VBF) 
processes qq — > q'q'Hf, remain nearly unchanged rela- 
tive to the corresponding processes in the SM. Thus, the 
corresponding SM cross sections and associated uncer- 
tainties described previously are also used here. In the 
FHM, direct decays to fermions are forbidden; the decays 
to W + W~ , 77, ZZ, and Zj account for nearly the entire 
decay width. For the mass range under investigation the 
W + W~ decay mode has the largest branching fraction. 
The branching fraction B(H f 77) is greatly enhanced 
over B(Hsm 77) for all m#, and the clean signa- 
ture and excellent mass resolution of this channel provide 



TABLE VIII: Observed and median expected (for the 
background-only hypothesis) 95% C.L. upper limits on the 
cross section times branching ratio cr(gg — > H) x B(H — >■ 
VF + VF~) from the combination of CDF and DO's Higgs bo- 
son search channels focusing on this production and decay 
mode, obtained using the Bayesian and CL B methods. 
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most of the search sensitivity for mu s < 120 GeV/c 2 . 
The analyses combined here seek Higgs boson decays to 
W + W~ , 77, and ZZ. Previous searches for a fermiopho- 
bic Higgs boson at the Tevatron excluded signals with 
masses smaller than 119 GeV/c 2 [88[; the expected ex- 
clusion was also mn f < 119 GeV/c 2 . The ATLAS and 
CMS Collaborations excluded mij, in the ranges 110.0— 
118.0 GeV/c 2 and 119.5-121.0 GeV/c 2 using diphoton 
final states [H] and in the range 110-194 GeV/c 2 by 
combining multiple final states |90j . 

The SM H — > 77 analyses are reoptimized as the 
kinematic distributions of the Higgs bosons, their de- 
cay products, and the particles produced in association 
with the Higgs bosons differ between the FHM and the 
SM. Events contain either an associated W or Z bo- 
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FIG. 24: (color online). Observed and median expected 
(for the background-only hypothesis) 95% C.L. Bayesian up- 
per production limits expressed as multiples of the SM4 (low- 
mass scenario) cross section as a function of Higgs boson mass 
from the combination of CDF and DO's Higgs boson searches 
focusing on the gg — > H production and H — > W + W~ decay 
modes. Uncertainties associated with theoretical cross sec- 
tion and branching ratio predictions are incorporated in the 
limit. The dark and light-shaded bands indicate, respectively, 
the one and two s.d. probability regions in which the limits 
are expected to fluctuate in the absence of signal. The red 
line shows the prediction for the signal rate in the high-mass 
scenario, divided by that of the low-mass scenario. 



son, or recoiling quark jets in the case of VBF and thus 
the transverse momentum (jpr) of the Higgs boson is 
on average greater than it is in the SM. The analyses 
combined here update previous searches for the Higgs 
boson in the FHM [HI H3- Similarly, SM searches in 
H — s- W + W~ channels cannot be interpreted directly 
in the FHM due to the different mixture of production 
modes. Signal contributions from gg — > Hf produc- 
tion to the MVA discriminant distributions are ignored, 
and the remaining contributions from other production 
mechanisms are scaled by the ratio of branching ratio 
predictions B(H f -> VV)/B(H SM -> VV). The exist- 
ing subdivision of channels based on the number of re- 
constructed jets accompanying the leptons and missing 
transverse energy in the event naturally optimizes the 
search within the FHM interpretation. Hence, the devel- 
opment of a separate set of analysis channels as in the 
case of Hf — > 77 is not required, though the MVAs are 
retrained. 

The combined limits on Higgs boson production nor- 
malized to FHM predictions obtained from both the 
Bayesian and CL S methods are listed in Table IIXI as a 
function of Higgs boson mass. The expected limits as- 
sume no Higgs boson production. The limits obtained 
using the Bayesian method are shown in Fig. I25I Fermio- 
phobic Higgs bosons in the mass range 100-116 GeV/c 2 
are excluded at the 95% C.L.; the expected excluded 
mass range is 100-135 GeV/c 2 . 



TABLE IX: Ratios of observed and median expected (for the 
background-only hypothesis) 95% C.L. upper limits on the 
production rate of a Fermiophobic Higgs boson relative to 
the FHM prediction as a function of the Higgs boson mass 
for the combination of CDF and DO's searches, obtained using 
the Bayesian and CL S methods. 
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FIG. 25: (color online). Observed and median expected (for 
the background-only hypothesis) 95% C.L. Bayesian upper 
production limits expressed as multiples of the FHM cross 
section as a function of Higgs boson mass for the combined 
CDF and DO searches. The dark and light-shaded bands in- 
dicate, respectively, the one and two s.d. probability regions 
in which the limits are expected to fluctuate in the absence 
of signal. 
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X. CONCLUSIONS 

The search for the standard model Higgs boson at the 
Tevatron is challenging due to the small expected signal 
and the need to accurately model large background con- 
tributions. We have developed advanced tools to search 
for the Higgs boson in the leading production and decay 
modes predicted by the SM and control the impact of sys- 
tematic uncertainties using constraints from the observed 
data. We have combined searches by the CDF and DO 
Collaborations for the standard model Higgs boson in the 
mass range 90-200 GeV/c 2 using Tevatron pp collision 
data corresponding to up to 10 fb" 1 of integrated lumi- 
nosity collected at y/s = 1.96 TeV. The results of searches 
focusing on the H -> bb, H -> W+W~ , H -> ZZ, 
H — > t + t~, and H — > 77 decay modes are included 
in the combination. The results are also interpreted in 
fermiophobic and fourth generation models. Fermiopho- 
bic Higgs bosons in the mass range 100-116 GeV/c 2 are 
excluded at the 95% C.L., and a SM-like Higgs boson in 
the presence of a fourth sequential generation of fcrmions 
is excluded in the mass range 121-225 GeV/c 2 at the 95% 
C.L. The SM Higgs boson is excluded, at the 95% C.L., 
from 90 to 109 GeV/c 2 , and from 149 to 182 GeV/c 2 . 
The expected exclusion regions in the absence of signal 
are 90-120_GeV/c 2 and 140-184 GeV/c 2 . The results of 
the H — > bb searches were validated through a measure- 
ment of the diboson (WZ + ZZ) production cross sec- 
tion using the same data samples and analysis techniques, 
treating those diboson processes as signal. The resulting 
diboson cross-section measurement is in agreement with 
the SM prediction. We observe a significant excess of 
events in the mass range between 115 and 140 GeV/c 2 . 



The local significance is maximum at mjj = 125 GeV/c 2 
and corresponds to 3.1 standard deviations, with a best- 
fit signal strength of 1.44j. ; 56 times the SM expectation. 
We also separately combined searches focusing on the 
H -> bb, H -> W + W~, H -> t+t~, and H -> 77 decay 
modes. The observed best-fit signal strengths obtained 
from each of these combinations arc consistent with the 
expectations for a SM Higgs boson at run = 125 GeV/c 2 . 
We performed tests of the compatibility of the observed 
excess with the expectations for the couplings of a SM 
Higgs boson and saw no significant deviations. 
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